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Abstract 

A frequency counter measures the input frequency V averaged over a 
suitable time r, versus the reference clock. High resolution is achieved 
by interpolating the clock signal. Further increased resolution is obtained 
by averaging multiple frequency measurements highly overlapped. In the 
presence of additive white noise or white phase noise, the square uncer- 
tainty improves from al oc 1/r 2 to &t oc 1/r 3 . Surprisingly, when a file of 
contiguous data is fed into the formula of the two-sample (Allan) variance 
a y( T ) = ^{^(Vk+i ~Vk) 2 } °f the fractional frequency fluctuation y, the 
result is the modified Allan variance mod o^(t). But if a sufficient number 
of contiguous measures are averaged in order to get a longer r and the 
data are fed into the same formula, the results is the (non-modified) Allan 
variance. Of course interpretation mistakes are around the corner if the 
counter internal process is not well understood. 
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1 Background 

Let v(t) = sin27w£ = sin[27Woo£ + the input signal, where tp(t) is the phase 
fluctuation, z^oo = l/2oo is the nominal frequency (the double subscript '00', as 
in foo, is used to avoid confusion with the 0-th term of a time series; thus, vqq 
is the same as vq commonly used in the literature, etc.), and u(t) = vqq + 
is the instantaneous frequency. Finally, let x(t) = 2 -if Voo the phase time fluc- 
tuation, i.e., the time jitter, and y(t) = 2 ^ 00 = V ~va°° ^ e fractional frequency 
fluctuation. The notation used in this article is the same of the general refer- 
ences about frequency stability and noise |BCC+72l lRut781 ICCI90I |Vig99| . 

Denoting with E{-} the expectation, the (classical) variance of y is <jy — 
E{[y — E{i/}] 2 }. In the presence of slow random phenomena, the variance de- 
pends on the measurement time and on the number of samples. This is related 
to the fact that the algorithm used is a filter in the frequency domain, whose 
lower cutoff frequency is set by the number of samples. Other variances are 
to be used, based on the idea that the estimator, clearly specified, has a lower 
cutoff frequency that blocks the dc and long-term components of noise. There 
results a variance which is a function of the measurement time r. Table^shows 
the spectral properties of the Allan variance and of the modified Allan variance, 
defined underneath. 



1.1 Allan variance (AVAR) 

Originally the Allan variance was introduced as a measurement tool for the 
frequency fluctuation of atomic clocks |A1166j . Given a stream of contiguous 
data y k averaged on a time r, the simplest variance is the (classical) variance 
evaluated on two samples, ct^(t) = ^[yk+i — yk] 2 - The estimated variance is 



Vk+1 Vk 



AVAR 



(1) 



and, expanding the time average 

r-(fe+2)T 



- y(t)dt-- 

T J{k+1)T T JkT 



(fc+l)r n2' 

y(t) dt 



(2) 



The above can be rewritten as 



a 2 y (r) =E 



+oo 

y{t) w A (t)dt 



(3) 
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Tabic 1: Noise types, power spectral densities, and Allan variances. 



noise 
type 


s v (f) 


Sy(f) 


S^p <— * Sy 




mod o 2 (r) 


white 
PM 


/),, 


hn f 2 


hn — H- 

"2 — -JI 


3f H h 2 -2 
(2tt) 2 
27tt/ h >1 


3f H T h 2 _-3 


flicker 
PM 


b-ir 1 


hif 


h 1 = b -^ 


[1.038 + 31n(27r/ ff r) 
X h < T~ 2 

X (27T) 2 T 


] 0.084 h lT - 2 
n>l 


white 
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b-2.r 2 
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h = b -£ 


\hor- 1 


jho r- 1 


flicker 
FM 


b- 3 f- 3 


h-if- 1 


h-i = 


21n(2) 


f ln(2) fc_ x 


random 

walk 

FM 


b-J- 4 


h-2.r 2 


h- 2 = ^ 


6 h - 2T 


0.824 (2 ; )2 /i- 2 r 


frequency drift y = 


Dy 




1 D 2 T 2 

2 


1 £)2 2 

2 -^ y 


Here vqq is replaced with vq for consistency with the general literature. 
Jh is the high cutoff frequency, needed for the noise power to be finite. 
The columns <j 2 (t) and mod o^(t) arc from Kro83. p. 79] (adapted). 



-^7 °<^<- 
-k r<t<2r 
elsewhere 



(Fig.© 



(4) 



which is similar to a wavelet variance. The weight function differs from the 
Haar wavelet in that it is normalized for power instead of energy. In fact, the 
energy of wa is 



£{wa] 



w 2 A (t)dt = 



1 



(5) 



while the energy of a wavelet is £ {•} = 1. 

In the frequency domain, the AVAR is similar to a half-octave bandpass 
filter with the peak at the frequency of o - . 



1.2 Modified Allan variance (MVAR) 

Another type of variance commonly used in time and frequency metrology is 
the modified Allan variance mod a 2 (r) |AB81I ILA84| . This variance was first 
introduced in the domain of optics Sny80| because it divides white phase noise 
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Figure 1: Weight functions for the two-sample variance and for the modified 
Allan variance. 



from flicker phase noise, which the AVAR does not. This is often useful in 
fast measurements. MVAR is also related to the sampling theorem and to the 
aliasing phenomenon jVZL98bl IVZL98a] because the trigger samples the input 
process at a rate 1/tq (see below). The MVAR is defined as 



mod cr^(r) 



e a 



-E -/ y(t)dt-- 



(i+n)T 



y(t) dt 
MVAR 



(6) 



with r = htq . 
The above is similar to a wavelet variance 

mod^(r) y(t) WM(t) dt 



(7) 



in which the weight function for tq <C t, or cquivalently for n 3> 1, can be 
written as 



W M = < ^ 



0<t<T 



V2t 2 

y(2t-3) r<t<2r 

-7I^(*- 3 ) 2r<t<3r 

elsewhere 



(Fig. IB 



(8) 



Once again, the weight function differs from a wavelet in that it is normalized 
for power instead of energy 



f + °° 1 
£{wai} = / w 2 M {t)dt= — 
J-oo 2r 



(9) 
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Figure 2: Rectangular averaging mechanism in simple frequency counters. 



Interestingly, it holds that £{wm} = ^£{wa}- This is related to the fact that 
the AVAR response to white frequency noise S y (f) = ha is cr^(r) = while 
the response of MVAR to the same noise is mod (r) = ^ . 



2 High-resolution frequency counters 

In this section the phase noise <p(t), and therefore the frequency noise y{t), is 
provisionally assumed to be a zero-mean stationary process. This hypothesis 
will no longer necessary for the direct measurement of the Allan variance. 

2.1 Classical reciprocal counters 

Traditionally, the uniform average over a suitable time interval r is used as an 
estimator of the frequency v. The expectation of v is therefore 

II estimator (10) 
(11) 

normalization . (12) 

Inside, the counter measures the time interval r = t^ — to between two zero 
crossings of v(t) spaced by N periods (Fig. Thus E{z/} = V = N/t. The 
averaging time tjj selected by the user is rounded to r = NT > tjj by stopping 
the measurement at the first zero crossing after that iy has elapsed. A variety 
of interpolation techniques |Kal04| can be exploited to avoid of the uncertainty 
-i- that results from the bare count of clock pulses at the frequency v c . The 
measurement of t is affected by the error xq — xn that results from the trigger 
noise and from the clock interpolator. Here, the reference clock is assumed ideal. 
Thus it holds that y = XN ~ Xa , With a state-of-the-art counter, the resolution of 
the interpolator in a single-event time interval measurement can be of 10~ n s. 



dt 



u(t)w n (t) 

-OO 

w u (t) = 



1/t < t < T 
elsewhere 



/+oo 
w n (t) dt = l 
-OO 
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Figure 3: Triangular averaging mechanism, implemented in some high-resolution 
frequency counters. 



Let us assume that xq and xn are independent and have identical statistical 
properties, and denote with a\ the variance of each. Under this assumption, 
which will be justified afterwards under more stringent conditions, the variance 
of r is 2a\ . Accordingly, the variance of the fractional frequency fluctuation is 

2 _ classical 
a v — T 2 variance ' ' 

The law <jy oc 1/t 2 is a property of the II estimator, i.e., of the uniform average, 
in the presence of white phase noise. 

2.2 Enhanced-resolution counters 

More sophisticated counters make use of the A estimator (Fig.[3J), which consists 
of a triangular-weight average. The counter takes a series of n measures T>i = 
N/ri delayed by iro = iDT, where ti = t,N+iD — Ud, i £ {0, ■ • • ,n — 1} is the 
time interval measured from the (i£>)-th to the (N + iD)-th zero crossings. The 
expectation of v is evaluated as the average 

71— 1 

E M = - ^ v i where V, = N/n . (14) 

i=0 



G 



Eq. (|14fl can be written as an integral similar to (|10fl . but for the weight function 
wu replaced with wa 

/+oo 
v(t)wA{t)dt A estimator . (15) 

-oo 

For To <C t, w\ approaches the triangular-shape function 

V < t < T 

w A (t) = { 2 - t/r r<t<2r (16) 
elsewhere 

w\(t)dt = l normalization. (17) 

-oo 

Nonetheless the integral l|15|) is evaluated as the sum <|14[) because the time 
measurements take place at the zero crossings. The measures F, are independent 
because the timing errors x^, k S {0, • ■ • , n — 1} are independent, as explained 
underneath. The counter noise is due to the interpolator noise, and to the noise 
of the input trigger. The samples of the interpolator jitter are independent 
because the interpolator is restarted every time it is used. The trigger noise 
spans from dc to the trigger bandwidth B, which is at least the maximum 
switching frequency of the counter. With modern instruments, B is hardly lower 
than 100 MHz, hence white noise is dominant. The autocorrelation function of 
the trigger noise is a sharp pulse of duration Tr ~ l/B. On the other hand, 
the delay tq is lower-bounded by the period Too of the input signal and by the 
conversion time of the interpolator. The latter may take a few microseconds. 
Hence in practice it holds that To 3> T^, and therefore the timing errors are 
independent. Accordingly, the variance of the fractional frequency fluctuation 
is 

2 _ 1 2o x classical /jo\ 
v n t-2 variance ' ^ ' 



At low input frequency, there is no reason for the delay r between overlapped 
measures Fj and Fj+i to be longer than Too, i-e., one period. Thus D = 1, 
to = Too, and n = N = vqqt. Hence Eq. ifT%)l is rewritten as 

2 _ 1 classical 

y ~ v 0Q r 3 variance ' ^ ' 

At high input frequency, the minimum delay tq is set by the conversion time 
of the interpolator. Hence the measurement rate is limited to vi measures per 
second, the number n of overlapped measures is n = vjt < voqt, and Eq. I|18f) 
turns into 

2 = 1M classical ( } 

y Vl r 3 variance v ; 

The law a 1 cx 1/t 3 , either (|19(l or J5DJ), is a property of the A estimator in the 
presence of white noise. This property is the main reason for having introduced 
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the A estimator in frequency counters. Yet it is to be made clear that the 
enhanced resolution is achieved by averaging on multiple measurements, even 
though overlapped, and that the measurement of a single event, like a start-stop 
time interval, can not be improved in this way. 

2.3 Understanding technical information 

Searching through the instruction manual of frequency counters and through 
the manufacturer web sites, one observes that the problem of the estimation is 
generally not addressed. When the counter is of the II type, the measurement 
mechanism is often explained with a figure similar to Fig. [JJ On the other 
hand, the explanation for the overlapped measurements in A-type counters is 
not found in the technical documentation. As a further clement of confusion, 
both counters provide one value every r seconds when programmed to measure 
over the time r. This can lead the experimentalist to erroneously assume that 
the estimation is always of the II type. 

The internal estimation mechanism can be understood from the formula 
for the "frequency error" , often given in the technical documentation. These 
formulae are of the form 

(H) ° V = \ ^W)li sge r + 2(^)? nter pol a tor (21) 

or 

(A) °y = ^2(St)l iggeI + 2((St)? ntor polator 

vqt i> w < vi 

n = < 

{ VlT V 0Q > VI 

where vj is of the order of 200 kHz. The actual formulae may differ in that 
uncertainty and noise of the reference frequency may be included or not; in that 
the factor 2 in the interpolator noise does not appear explicitely; and in other 
details. 

The terms inside the square root of (5T) and (55) come from independent 
white noise processes, as explained in Sections |2.1l Thus, one can match Eq. I|21[) 
to and Eq. (55) to (|18|) . Consequently, the presence of a term r in the 

denominator reveals that the counter is of the II type, while the presence of the 
term r^fn or Ty/r reveals that the counter is of the A type. 

2.4 Examples 

Two instruments have been selected with the sole criteria that the author is 
familiar with them, and that they are well suitable to show how to the II and 
the A estimators can be identified in the instruction manual. 



Stanford Research Systems SR-620 (n estimator). In the instruction 
manual |Sta04l p. 27], the RMS uncertainty is called resolution and given by the 
formula 



RMS 

resolution 
(in Hz) 



frequency 
gate time 



\ 



(25 ps) 2 



( short term \ x / gate A 
I stability J ytimc J 



2x 



trigger 
jitter 



N 



(23) 



The above <|23|) matches our notation with 



RMS resolution a v = 
frequency voo 
gate time r 



voo&y (classical variance) 



The numerator inside the square root is the square single-shot time deviation 
2er 2 . It includes the inherent resolution of the counter, 25 ps, ascribed to the 
interpolator; the phase noise of the frequency standard; and the equivalent input 
noise (signal plus trigger) divided by the input slew rate. The denominator 
"N" is the number of measurement averaged, not to be mistaken for N of our 
notation. If the measurement time is equal to the gate time, it holds that N = 1. 
Eq. divided by Uqq and squared is of the same form of Eq. (|T3|l . thus of Eq. 
<|21|) . with dy oc . This indicates that the counter is of the II type. 



Agilent Technologies 53132A (A estimator). The instruction manual 
|Agi98[ pp. 3-5 to 3-8] reports the RMS resolution given by the formula 



RMS 
resolution 



f frequency \ ^ 
yor period I 



2 X (trig 



;ger error 



t 



jitter 



(gate time) X no. of samples gate time 



i — 225 ps 



t 

i jitter 



3 ps 

number of samples 



(24) 



(gate time) x (frequency) 
(gate time) x 2xl0 5 



for / < 200 kHz 
for / > 200 kHz 



The above l|24(l matches our notation with 

v m a y or a T - 



RMS resolution 
frequency 
period 
gate time 
no. of samples 

/ 



OV = 
Too 

T 

n 



Too&y 



(classical variance) 



The term a^Tmc of (EU is due to some internal phenomena not analyzed here, 
for we provisionally neglect it. The numerator inside the square brackets is the 
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single-shot time deviation 2tr^. After squaring and normalizing, Eq. I|24|) is of 
the same form of (I18|) . thus of (|22|l . The counter switches averaging mode at 
the input frequency of vi — 200 kHz. In fact, it holds that 




t'ooT t'oo < 200 kHz 

rx2x!0 5 t/ 00 > 200 kHz 



Accordingly, Eq. {2U turns into {Hi f° r v < 2xl0 5 Hz, and into (PJJJ beyond 
2xl0 5 Hz. The variance is oc in both cases, for the counter is of the A 
type. 



Remarks. The resolution enhancement of the A average can only be used 
with stationary phenomena, not with single events. This is easily seen in the 
table underneath. This table compares the two counters, SR-620 and 53132A, 
under the simplified assumptions that the reference clock is ideal, and that the 
high slew rate makes the trigger noise negligible. 

counter single-event start-stop frequency vqq = 100 kHz 
type time interval TI = 1 s averaged on r = 1 s 

SR-620 o- = 2.5xl0" n a y = 2.5X10" 11 

53132A (T = 9xl0- 10 a y = 5.8xl0~ 12 

The superior resolution of the SR620 interpolator provides higher resolution for 
the measurement of a single event in start-stop mode. On the other hand, the 
highest resolution is obtained with the 53132A in the measurement of stationary 
signals, due to the A estimator. 



3 Frequency stability measurement 

Let us first observe that it holds 

w A (t) = -=\w n {t - t) - w u {t)] (Fig.© (25) 



V2 
1 

V2 



w M (t) = ^= to A (t- t) - w A (t) (Fig.© . (26) 



3.1 Il-type counters 

AVAR. Let us get a stream of data yt from the output of a n-estimator 
counter, measured over a base time slot tb with zero dead time. The superscript 
"(i)" rc fcrs to the averaging time 1 x tb- Feeding this stream into Eq. (JJJ, one 
gets Eq. J2J) with r = tb- Then, after averaging contiguous data in groups of 
to, one gets a smaller file of averaged over the time ttltb- Feeding this 

new file into Eq. QJ, one gets Eq. J5J) with r = jtitb- This is exactly what one 
expects. It is a common practice to plot of u^(r) in this way, using a single data 
stream and m in powers of 2. 
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AVAR 
w n (t) 

W n (t-T) 
W A (t) 

x 2x 



MVAR 



w A (t) 


\ • 1/T 

^ time 


W A (t-T) 


\ 1/T 1 

: - 


w m(0 


\ 1 
•• +1/(V2t) 




; l/( . 2 T) l' 



t 2t 3t 

Figure 4: Relationships between the weight functions. 

MVAR. The data file y^ can also be fed in Eq. JHJ. In this case the mea- 
surement time tb of the counter is the delay tq of ©■ The variance is evaluated 
at t = titb set by the choice of n. Yet, it is desirable that n> 1. 



1/T 



time 



1/T 



+1/(V2t) 



-l/fV2Tl 



3.2 A-type counters. 

The attention of the reader is now to be drawn to a subtilty in the use of a 
A-estimator to measure the Allan variance. While the counter provides one 
value of v every tb seconds, two contiguous windows w^{t) and w\(t — tb) are 
overlapped by tb- That is, the falling side of w\(t) overlaps the rising side of 
w\(t — tb)- Thus Vk+i — Vk, hence y^+i ~Vk^ i s the frequency averaged with 
the wm window. The practical consequence is that, feeding the file of such 
into Eq. (JTJ (AVAR), one gets exactly Eq. © with t = t B - That is, MVAR 
instead of AVAR. 

A longer measurement time r is obtained by averaging contiguous data in 
groups of to. This process yields a smaller file of y k m \ averaged over tots. Yet 
the measurements are overlapped, for the weight function, which is the triangle 
w\ for to = 1, turns into an isosceles trapezium for to > 1, and asymptotically 
into the rectangle wjj for large to. A graphical proof is proof is given in Fig. [5J 
The practical consequence is that, feeding the file of |k into Eq. 0J, one gets 
an "odd" variance for with small m > 1, i.e., a variance that is neither AVAR 
or MVAR; and again the AVAR for large m. 

Unfortunately, the conversion between II and A estimates is impossible with- 
out combining a large number of data. It is therefore impossible to measure the 
AVAR with a A-type counter without increasing r, that is, r = vtltb with a 
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Figure 5: The average of m contiguous A measures converges to a II estimate. 

sufficiently large value of m. 

Finally, we focus the attention to two asymptotic classes of measurement, 
namely short-term stability and long-term stability, analyzed underneath. 

Long-term stability. These measurements are characterized by a large value 
of m. Hence the averaging function converges to wu, thus the data file can 
be fed into the AVAR formula JJJ without mistakes. This is a fortunate outcome 
for two reasons, namely 

1. AVAR is preferable to MVAR because, for a given duration of the time 
record, it provides a more accurate estimate at large r, 

2. the rejection of white phase noise of the MVAR is not necessary in long 
term measurements. 

Short-term stability. The experimentalist interested in short-term stability 
appreciates the rejection of white noise of the MVAR, while the longer duration 
of the experiment, as compared to the measurement of the AVAR at the same 
r is not disturbing. In this case, the bare mean can not be used to combine 
contiguous values in order to get t = ttitb- The values must be weighted 
proportionally to the triangular staircase sequence {l,2,...,[m/2],...,2,l}, 
so that the equivalent weight function is an isosceles triangle of width It. A 
file of such measures fed into gives the MVAR evaluated at r = ttitb- The 
nuisance is that this triangular-shape average is only possible for odd m. 

3.3 Other methods 

For long-term measurements, the total variance TotVar |GHP99) is progressively 
being used as an estimator of the AVAR. The importance of the TotVar resides 
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in that, given a record if measured data, it exhibits higher accuracy at larger r. 
Yet TotVar is based on rectangular averages. Consequently, the raw readout of 
a A estimator can not be fed in the formula for TotVar without interpretation 
mistakes. 

Another useful method, often called picket-fence, consists of the absolute 
timing of the zero crossings versus a free-running time scale. In this case, there 
is no preprocessing inside the counter, for there can be no ambiguity in the result 
interpretation. Absolute-time data can be used to calculate AVAR, MVAR, 
TotVar an d other variances. The concept of picket-fence was first proposed in 
|BCC+72I Eq.(79)]. The n, it was studied extensively as an independent method 
for the measurement of AVAR and MVAR |Gre89l lGre97] . 
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